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ABSTRACT 

The recent discovery of hidden non-axysimmetric and disk-like structures in bright 
Virgo dwarf elliptical and lenticular galaxies (dE/dSph/dSO) indicates that they may 
have late-type progenitors. Using N-body simulations we follow the evolution of disk 
galaxies within a ACDM cluster simulated with 10 7 particles, where the hierarchical 
growth and galaxy harassment are modeled self-consistently. Most of the galaxies 
undergo significant morphological transformation, even at the outskirts of the cluster, 
and move through the Hubble sequence from late type disks to dwarf spheroidals. 
None of the disks are completely destroyed therefore they can not be the progenitors 
of ultra compact dwarf galaxies (UCDs). The time evolution of the simulated galaxies 
is compared with unsharp-masked images obtained from VLT data and the projected 
kinematics of our models with the latest high resolution spectroscopic studies from 
the Keck and Palomar telescopes. 

Key words: methods: N-body simulations - galaxies: clusters: general - galaxies: 
dwarf - galaxies: evolution 



1 INTRODUCTION 

Early type dwarfs (dE and dSO) are the most common 
type of galaxies in the nearby universe, yet their origin and 
nature are still unknown. These galaxies are found nearly 
exclusively near bright galaxies and in groups or clusters. 
Their continuation or separatio n from the brighter elliptical 
sequence is currently deba ted dGraham. Jerien fc Guzman! 
2003; iGavazzi et all l2004h . Formation scenarios include 
mainly two different hypotheses: 1) dEs and dSOs are 
primordial galaxies, 2) they are the result of a morpho- 
logical transformation of spiral and irregular galaxies 
accreting into the cluster. This latter idea is supported 
by some observations, i.e the relative number density of 
dwarf e arly type galaxies increase s with the local galaxy 
density iFerguson fc Sandagejl99rh . suggesting that the en- 
vironm ent drives galaxy evolution. iKormendv fc Freeman! 
( 2004) investigated the systematic properties of dark 
matter halos in late-type and dwarf spheroidal galaxies 
concluding that they form a single p hysical sequence as a 
functi on of the dark matter core mass. iMoore. Lake fc Katzl 
(1998) have shown that galaxy harassment in clusters 
can transform spirals into spheroidals. Harassment can 



explain the morphological evolution of the small spiral and 
irregular galaxies ob served in clusters at redshift z ~ 0.4 
toressler et al.| | l994h into dwarf ellipticals in the nearby 
universe" ilMooreeT^l]ll996l) . Moreover the radial velocities 
of early type galaxies in the Virgo cluster seem to indicate 
that these galaxies are not an old cluster population 
JConselice. Gallagher fc Wvseil20oJl but originate from in- 
falling field galaxies. According to this evolutionary scenario 
we expect to find nearby cluster galaxies that are currently 
undergoing morphological transformation and retain part 
of their disk nature. Recent observations confirm that early 
type galaxies have a bro ad range of photometric and kine - 
mati cal characteristics. [Jerjen. Kaln^s fc Binggelil (l200Ct) 
and iBarazza. Binggeli~fcJferlenn2002ri discovered hidden 
spiral structures and bar features in five bright dEs in the 
Virgo cluster, concluding that the fraction of early type 
dwarfs hosting a disk component could be larger than 20%. 
Spirals and disks were al so observed within dwarf spheroida l 
galaxies in the Coma j Graham 1 Jeri en fc Guzmanl 120031) 
and Fornax clusters foe Riicke et alJ l200lt) ~ Furthermore, 
the degree of rotational support is found to vary from zero 
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12001 120031: Ivan Zee. Skillman fc Havnesll2004li . 

The aim of this work is to follow the evolution of disk 
galaxies orbiting in a cluster environment using high res- 
olution N-body simulations and to compare the final ha- 
rassed remnants with the latest photometric and spectro- 
scopic data. Previous studies have suffered from low reso- 
lution and they used idealised cluster models for the initial 
conditions. We aim to use sufficient resolution that we can 
follow the detailed morphological evolution and we evolve 
the galaxy models within a cluster selected from a cosmo- 
logical simulation. In addition, the model has cuspy CDM 
halos and structural parameters as expected in the concor- 
dance cosmolo gical model. This study complements that of 
iGnedinl i2003T) who studied the evolution of luminous early 
type disk galaxies within a cosmological context. 

The paper is organised as follows. In Section 2 we 
present the main characteristics of the galaxy models and 
the ACDM cluster used. The results of the three dimensional 
simulations are illustrated in Section 3, while in Section 4 
and 5 we report a projected photometric and kinematical 
analysis of the simulated remnants and the comparison with 
observations. 
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Figure 1. The rotation curve of the initial galaxy model. The 
stellar and dark matter contributions to the total rotation curve 
arc indicated. 



2 SIMULATIONS 

All the simulations have been carried out using PKDGRAV, 
a parallel multi stepping N-b ody tree-code designed for high 
force accuracy JStadeH200lD . 

The galaxy model is a multi-component system with a 
stellar disk embedded in a spherical dark matte r halo and 
was co nstructed using the technique described bv lHernauisd 
(1993). The dark matter distribution init ially follows a 
NFW iNavarro. Frenk fc White! Il99rj . Il997h profile, which 
is adia batically contrac ted in response to baryonic infall 
iSpringel fc WhitellT9 99') . The disk has an exponential sur- 
face density profile of the form: 

Z(R) = ^h^p(-R/R d ), (1) 

where Md and Rd are, respectively, the disk mass and radial 
scale length (in cylindrical coordinates), while the thin ver- 
tical structure is characterised by the scale height Zd which 
sets the "temperature" of the disk 

p d (R,z) = ^sech 2 ( Z /^). (2) 

^Zd 

The structural parameters of the disk and the halo are cho- 
sen so that the resulting rotation curve resembles that of a 
typical bulgeless late-type (S c/Sd) disk galaxy dCourteaul 
Il997l: iPersic fc Saluccl Il997fl . The model pa rameters are 
initial ised following the same procedure as iMaver et alJ 
( 2002 ), which was based on the galaxy formation model of 
iMo. Mao, fc White! (Il99ct) . The mass within the virial radius 
was set equal to 7 x 10 10 Mq and the fraction of mass in the 
disk is ~ 6%. The contribution of the different components 
to the global rotation curve, assuming a disk scale length 
Rd = 1.5 kpc and a concentration c = 10 (where c is defined 
as c = r V ir/r 3 , with r„j r and r 3 respectively virial and scale 
radius of the NFW halo) is plotted in Fig. 1. The halo spin 
parameter, which sets the disk scale length in our model- 
ing, is A = 0.045, where A relates the angular momentum J 




Figure 2. Q and X2 parameters as function of radius (in units 
of the disk scale length Rd)- 



and the total energy E of a system with virial mass M v i r 
through the relation A = J\E\ 1/2 G- X M'^J 2 . The model has 
a central surface brightness /j,b = 22magarcsec~ 2 , assum- 
ing an initial V-band mass to light ratio ~ 4 GEha et al 
2002 and & B — V color of 0.77 Ivan Zee. Barton fc Skillma: 
2004). The stellar disk satisfies the Toomre llToomrelll9~ 
stability criterion, which requires, for a rotational supported 
disk 



Qstar{R) — 



<?nk 
3.36GXL. 



> 1, 



(3) 



where or is the radial velocity dispersion, k is the local 
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Figure 3. Unsharp masking technique applied to the stellar component of a simulated face-on galaxy (GAL5 of Table 1). The box size 
is ~ 20 kpc. Left frame: smoothed stellar surface density produced from the the simulation output using BEAM by J. Stadel. Middle 
frame: the result of the unsharp masking on the former image, which reveals presence of a bar and a spiral pattern. Right frame: a noise 
filter is previously applied to the original surface density map, increasing the intensity of the gray channel for random selected pixels. As 
a result in the final unsharp image the low density structures disappear and only the bar and a ring feature are still visible. 




Figure 4. Morphological transformation (face-on projection) of a disk galaxy into a dwarf elliptical (GAL8 of Table 1). For each frame 
the redshift is indicated. The box size is ~ 20 kpc. The step like structures visible at the edges of the galaxies are due to features present 
in the low density regions of the original smoothed stellar surface density maps, which are amplified by the unsharp masking. 



epicyclic frequency and E 3 the stellar surface density. The 
efficiency of swing amplification of a disk perturbation with 
m-fold symmetry is governed by a co mbination of Q and 
the parameter X m = k 2 Rj '(47rmG£ s ) llToomrelll98lh . The 
radial profiles of Q and X2 are indicated in Fig. 2. 

Each galaxy is modeled with 10 6 particles, 10 s of which 
are in the disk. The gravitational softening is set to 0.5 



kpc for the dark halo and 0.1 kpc for the disk particles. 
The galaxies were evolved for ~ 8.4 Gyrs in a ACDM clus- 
ter with a mass at z = comparable to that of Virgo 
{M v i r = 3.1 x 10 14 M©) and a virial radius R V i r of 1.8 
Mpc. The cluster, selected from a 300 Mpc cube simu- 
lation JPiemand. Moore fc Stadel l2004bT) with cosmologi- 
cal parameters Qa = 0.732, Q m = 0.268 and as = 0.9 
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redshift z are expressed in terms of the Hubble constant 
H(z) through the relation 



Figure 5. Evolution of a representative galaxy (GAL7 of Table 
1) from redshift 0.5 to the present time. Top panel: strength of 
the bar. Bottom: buckling amplitude. The bar forms at z ~ 0.33 
and strongly buckles at z ~ 0.2. 




Figure 6. Evolution of the stellar surface density of GAL7 from 
z = 0.3 (left panels) to z ~ (right panels). The color scale is 
logarithmic. The top panels represent the face on projection while 
the corresponding edge on projection is shown on the bottom. The 
buckling phase, which has a peak at z = 0.2, does not destroy the 
bar but makes the central vertical structure of the galaxy thicker. 



flSpergel et all20 03). was resimulated at high resolution and 
corresponds to the intermediate resolution cluster D9 of 
iDiemand. Moore fe Stadel <2004a^ . The number of particles 
within the virial radius is ~ 6 x 10 6 and the softening length 
is 2.4 kpc. 

At 2 = 0.5 we replaced 20 random cluster particles (10 
within the virial radius and 10 within 20% of the virial ra- 
dius) with the high resolution galaxy model. Only one galaxy 
model is used, in order to isolate the effects of orbit on the 
evolution. We assigned to the center of mass of each galaxy 
the same orbital position and velocity of the replaced parti- 
cle. We replace particles rather than halos so as not to bias 
the results by omitting "overmerged" halos. Internal veloci- 
ties and lengths were rescaled according to the redshift and 
cosmology. In particular the comoving velocities at a given 



H{z)x, 



(4) 



where r and x are the physical and comoving coordinate and 
a the scale factor 1/(1 + z), respectively. 



3 GALAXY EVOLUTION 

Among the initial sample of 20 galaxies, only 13 are still 
within the cluster virial radius at z = 0, while the oth- 
ers lie on bound orbits between Rvir and 1.5 Rvir- This 
result is consistent with the fact that about half of the 
halos which are presently at the outskirts of the clus- 
ter (up to 2R v j r ) have a peri center smaller than R V i r 
jMoore. Diemand fc Stadelll2004l) . 

The evolution of many of the galaxies within the clus- 
ter is quite violent and due to a sequence of strong gravi- 
tational encounters with substructures and with the global 
cluster potential. Most of the central galaxies lose a signifi- 
cant fraction of stars and undergo a complete morphological 
transformation from disks to spheroidal systems. In order 
to highlight the presence of hidden structures we applied to 
the stellar component of the simulation results the unsharp 
masking technique described in Barazza et al. 2002. This 
method consists of the following steps. The intensity of each 
pixel is replaced with the mean intensity of a given area 
around the pixel producing a smoothed image. The linear 
dimension of this area has to be close to the typical scale 
of the features that we want to uncover. The original image 
is then divided by the smoothed one. Fig. 3 illustrates the 
results of the unsharp masking applied to a galaxy of our 
sample, where the ESO's MIDAS task FILTER/SMOOTH 
is used for the image smoothing with a smoothing length of 
30-40 pixels, corresponding to ~ 2 kpc. 

Fig. 4 represents the different phases of the evolution 
of a galaxy transforming into a dE. The first stage of the 
evolution is characterised b y the formation of a strong bar 
jMoore. Lake fc Katdll998h and an open spiral pattern in 
the disk. The spiral arms are easily stripped while the re- 
maining material forms a ring structure around the bar. The 
bar+ring phase is quite stable and is apparent in the final 
state of several of the remnant galaxies. Often the ring is not 
directly visible in the simulation output, but is clearly seen 
after applying the unsharp masking. When gravitational 
heating due to tidal interactions removes the remaining spi- 
ral and ring features, only a naked bar remains. The bar usu- 
ally undergoes a strong buckling instability and as a result 
the central part of the galaxy becomes more spherical. The 
buckled bar is subjected to further tidal heating and loses 
mass from its edges, becoming rounder the more encounters 
it suffers and the longer it evolves near to the cluster centre. 
Fig. 5 illustrates the evolution of bar and buckling instabili- 
ties for a representative galaxy which is orbiting close to the 
center of the cluster at z = 0. Abar and Abuchling are the bar 
and buckling amplitudes, defined respectively as the Fourier 
decomposition of the phas e of particles in the xy pl ane and 
of their vertical positions dSparke fc Sellwoodi ri987) 
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Figure 7. Three different final states of the harassed galaxies (unsharp nosified images). From the left to the right: bar-like galaxy 
(GAL7 of Table 1), spherical spheroid (GAL6) and spiral galaxy (GAL3). The addition of noises often hides ring structures surrounding 
the central bar, which is the case of the first remnant on the left. In order to point out the characteristic strucure of the remnants we 
used two physical scales: the two images on the left have a box of ~ 12 kpc while the spiral galaxy on the right has a diameter of ~ 20 
kpc. 




Figure 8. Tidal tails of GAL6 (Table 1) generated during a peri- 
centric passage (unsharp nosified image). These smooth features 
are visible only for a short time (typically about 0.2-0.3 Gyr) since 
they are heated by continued encounters with cluster substruc- 
tures. 
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Fig. 6 shows a face and edge on projection of the same galaxy 
during the buckling phase where the bar appeared at z = 
0.3. 

Depending on the orbit and on the number of close en- 



Figure 9. Intracluster light produced by tidal processes. The 
color contrast is chosen in such a way as to point out the low 
density regions (light green, corresponding to a surface brightness 
^<,30 mag arcsec - 2 ) and the traces of stellar streams. 



counters that the galaxy suffers, the morphological transfor- 
mation of the stellar component can be complete and pro- 
duce a pure spheroidal system, or can lead to the formation 
of a more elongated structure, which still retains disk and 
bar features (Fig. 7). For a few galaxies orbiting closest to 
the cluster center the tidal effects from the cluster poten- 
tial become dominant and the evolution is very rapid. Fig. 
8 shows the tidal tails generated by a strong impulsive en- 
counter. The stripped mass is a large fraction of the total 
stellar mass and the surface density in the tails is compara- 
ble with that of the outer disk. We do not expect that this 
kind of structures will be easily observable in clusters, since 
their life time is quit short (« 0.2 Gyrs) due to numerous 
encounters with other substructures which heat and dissi- 
pate the stripped material. Stripped stars tend to create a 
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Figure 10. Final flattening of the remnants plotted versus the 
fraction of stellar mass within the tidal radius at z = 0. 
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Table 1. Final state of the remnants. The second column indi- 
cates, when available, the apocentric and pericentric distances of 
the last orbit. For each galaxy we list three dimensional (spher- 
ically averaged) tidal rt and effective r e radius (in kpc), dark 
matter M^ and stellar M s mass (in units of IO^Mq). The tidal 
radius is calculated using SKID (Stadcl 2001), while r e is the ra- 
dius containing half of the light. The last two columns indicate 
the axial ratios (respectively intermediate and short axis to major 
axis) measured within 2r e . 



diffuse distribution around the remnants (see Fig. 9) and 
in regions close to the center it is almost impossible to dis- 
tinguish well defined stellar streams, even from the debris 
of a small fraction of the true number of dwarf galaxies in 
clusters. 

Table 1 summarises the main properties of the remnant 



galaxies. The initial model has a virial radius of 85 kpc and 
the effective radius of the stellar distribution is 2.5 kpc, de- 
fined as the radius containing half of the light of the galaxy. 
While almost all the galaxies lose more than two thirds of 
the dark matter halo, galaxies located at outskirts of the 
cluster do not experience a significative amount of star loss, 
even if their stellar structure can be perturbed producing 
open spiral patterns and asymmetric features (last frame on 
the right of Fig. 7) . For these remnants the reduction of the 
effective three-dimensional radius r e is related to the for- 
mation of a bar that increases the central stellar density. 
On the other hand, galaxies orbiting closer to the cluster 
center lose up to 90% of their stellar mass and the decrease 
of r e corresponds to a real decrease in the size of the stel- 
lar component. In a few cases the galaxy loses much more 
dark matter than stars and becomes baryon dominated in 
the central region. This is due to the fact that the orbits 
of the dark matter particles in the halo are more eccentric 
than those of the star particles and therefore reach the tidal 
radius more easily. Table 1 shows that the loss of stellar 
mass is not a simple function of the mean orbital radius or 
apocentric/pericentric distance from the cluster center, con- 
firming the importance of harassment. The last two columns 
of Table 1 give the mean axial ratios within 2r e . There is 
a clear correlation between the c/a ratio (the flattening of 
the stellar component) and the final stellar mass. As shown 
in Fig. 10, massive stellar remnants have smaller c/a values 
and retain more of their initial disk nature, while galaxies 
that have lost most of their stars tend to be more spheroidal. 
The final state is always quite prolate, with the exception 
of GAL6 and GAL8 (the first two points in the upper left 
corner of Fig. 10) which are close to spherical. In the case 
of the less perturbed galaxies this is due to the fact that 
the radius 2r e defines an area roughly corresponding to the 
bar region, while most of the other galaxies retain part of 
the radial anisotro py that originated during the bar phase 
jMaver et al.ll200ll) . 



4 PHOTOMETRIC ANALYSIS 

In order to compare the results of our simulations with ob- 
servat ions and in particular with lBarazza. Binggeli fc Jerienl 
(2003), we selected the remnants within the virial radius at 
z = (first 13 galaxies of Table 1) and projected them along 
a random line of sight. 

The isodensity contours of the projected disks were fit- 
ted with ellipses using the MIDAS task FIT/ELL3. The el- 
lipses centers, ellipticity and position angle were considered 
as free parameters of the fitting procedure. The isophotal 
contours of a representative galaxy (GAL9) are shown in Fig 

11. Ellipticity and position angle profiles are plotted in Fig. 

12, where the independent variable is the equivalent radius 
r p = a p b p , and a p and b p are the major and minor axis 
of the projected remnant. Even though there is not a well 
defined trend in the ellipticity profiles, galaxies which pre- 
serve a bar component are generally characterised by large 
ellipticities in the central 2-3 kpc and by a drop in the outer 
regions. Often the inner 0.5 kpc exhibits a very low ellip- 
ticity corresponding to a small nuclear component produced 
by buckling instabilities, whereas the external drop is more 
likely due to disk heating. Rounder systems have a wider 
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Figure 12. Each set of three vertical panels represents the photometric characteristics of a galaxy from Table 1: B-band surface brightness 
(top panel) , ellipticity (center) and position angle profile (bottom) . The coordinate r p is the equivalent radius r p = \J a p b p , where a p and 
b p are the major and minor axis of the projected remnant. The Sersic model is superimposed on the surface brightness profile (red solid 
line) and the corresponding best fit parameters are given in the top right of the panel. The central surface brightness fio is expressed in 
magarcsec - 2 , rrj is in kpc and the shape parameter n is a pure number which represents the deviation of the profile from an exponential 
law. 
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Figure 12. continued. 
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Figure 12. continued. 
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Figure 12. continued. 




Figure 13. B-band surface brightness, ellipticity and position angle profiles of the face on and end on projections of GAL7. 
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Figure 11. Elliptical contours of GAL9. A strong isophotal twist 
is evident in the outer region (see also Fig. 12). 



Figure 14. Surface density distribution of the two remnants pre- 
senting a nearly face on spiral pattern (GAL3 on the left and 
GAL12 on the right). The physical size of both the images is 
~ 12 kpc. 
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Table 2. Effective radii (in kpc) of the projected remnants calcu- 
lated using the integrated light curves. In most of the cases they 
do not deviate considerably from the values obtained for the three 
dimensional galaxies. A significative exception is GAL13, a disk 
system viewed almost edge on. 

range of profiles, with eccentricities that may even increase 
with radius, as in the case of GALl. 

The next step of this analysis consists in calculating the 
surface brightness profiles of the projected remnants. In or- 
der to derive a luminosity we need to assign a mass to light 
ratio. We expect that the ISM would be rapidly stripped by 
ram pressure effects such that star formation is truncated 
shortly after the galaxies enter the cluster environment. In 
this case, roughly 8 Gyrs from the beginning of the sim- 
ulation, the stellar mass to li ght ratio will have increased 
by about a factor of 1.5 — 2 iMaver et alJl2001f) . For this 
reason we used a B-band M/L ~ 6 for the final remnants. 
The cumulative light profiles were determined by integrating 
the light of each galaxy in elliptical apertures with increas- 
ing equivalent radius (MIDAS task INTEGRATE /ELLIPS) 
and fixed center, ellipticity and position angle. For these 
parameters we adopted their mean values within a radius 
corresponding to half of the disk size, in order to avoid un- 
certainties in the isophote fitting procedure due to the low 
density external regions. In Table 2 we list the effective ra- 
dius of each projected galaxy defined by r ep which is the 
radius at which the cumulative light profile reaches half of 
the value that it has at the isophotal level of 28 mag arcsec -2 
(upper B-band magnitude limit in Barazza et al. 2003). 

The surface brightness profiles, obtained by differenti- 
ating the cumulative light curve with respect to the equiv- 
alent radius r v (Barazza et al. 2003), are plotted in Fig. 
12. The oscillation in the profiles of GAL3 and GAL12 re- 
veals the presence of spiral arms viewed almost face on (Fig. 



14), which enhance the local stellar density and produce 
abrupt changes in the position angle and eccentricity pro- 
files. The surface brightness profiles between 0.1 kpc (soft- 
ening length for the stellar component) and the isophotal 
level of 27 mag arcsec -2 were fitted with a Sersic profile 

<7(r) = aoe ( - r/r ° )<1/n \ (7) 

where a is the surface brightness and ero (central surface 
brightnes), ro (scale length) and n (shape parameter) are 
the free parameters of the fit. For n — 1 the Sersic model 
corresponds to an exponential law, while a shape parameter 
n = 4 is indicative of a de Vaucouleurs profile. Converting 
surface brightnesses to magnitudes, Equation |7| becomes 

H(r) =^o + 1.086(r/r ) (1/n) , (8) 

where is the magnitude corresponding to ctq. 

The fitting curves and parameters of those galaxies that 
it was possible to fit with a Sersic model are indicated in 
Fig. 12. The range of surface brightness parameters agrees 
quite well with the observations of dwarf galaxies in Virgo. 
The central surface brightness of our remnants is on av- 
erage slightly higher than the values observed in B-band 
by Barazza et al. 2003, who fitted Sersic profiles over simi- 
lar radial intervals, the scale lengths obtained from the fit- 
ting procedure are comparable with their results. Also note 
that the quoted central surface brightness values are Sercic 
fits extrapolated to r — 0. The actual values at our reso- 
lution limit are somewhat fainter (Table 3). Galaxies with 
the smallest values of the shape parameter n are not well 
fitted by a Sersic profile in the central part. This is the case 
of GAL4 - apparently a pure spheroidal system with high 
ellipticity, in reality characterised by a disk component seen 
edge on - and of GAL7, which is still an asymmetric bar like 
system with a large central nuclear component with small 
ellipticity values. For GAL7 we also show surface brightness 
profiles, ellipticity and position angles for the face on and 
the end on projections (Fig. 13). The smallest central ellip- 
ticity is associated with the end on projection, where the 
bar is viewed along the major axis. Only the less massive 
remnants GAL6 and GAL8 have surface brightness profiles 
approaching to exponentials (n ~ 1). 
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Figure 15. Projected kinematic profiles of the remnants within the cluster virial radius at z = 0. The line of sight velocity (left panels) 
and the velocity dispersion (right panels) are plotted as a function of the radial distance along the major axis. 
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Figure 16. Projected kinematic profiles of GAL6 and GAL8 for i = 0° (face-on orientation), i = 45°, % = 90° (edge-on) and i = 90° 
with the three dimensional major axis aligned with the line of sight (end-on projection). The line of sight velocity (left panels) and the 
velocity dispersion (right panels) are plotted as a function of the radial distance along the projected major axis. 
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Table 3. Parameters of the Sersic fit: central surface brightness 
(in mag arcsec -2 ), surface brightness at the softening radius, scale 
length (kpc) and shape parameter. 

5 KINEMATICS 

We attempt to measure the kinematics of the remnants in 
a similar way as previous authors measuring dwarf galaxies 
in the Virgo cluster. We measure the kinematics of the 13 
projected galaxies analyzed in Section 4 using a 0.3 x 12 
kpc s lit, which corresponds, at the distance of Virgo (15.3 
Mpc. lFreedman et al.l J200ll) 1 to 3" x 40". This radial range 
is comparable with the slit length adopted by van Zee et 
al. 2004 and allows us to observe in most of the galax- 
ies the turnover of the rotation curve, while the thickness 
of the slit was set to the softening length of 0.1 kpc. The 
slit was positioned along the major axis of the projected 
remnant centered on the galaxy and star particles within 
the slit area were binned into a grid. For each grid inter- 
val with more than 15 particles - excluding thus the most 
weakly bound particles at the edges of the system - we cal- 
culated the mean particle velocity (relative to the center of 
mass of the galaxy) and velocity dispersion along the line 
of sight. For the two smallest remnants GAL6 and GAL8 
we increased the number of bins in order to obtain more 
precise velocity curves. The resulting kinematic profiles are 
indicated in Fig. 15, where for each galaxy the line of sight 
velocity (left panels) and velocity dispersion (right panels) 
are plotted as a function of the radial distance along the 
major axis. In most of the cases the remnants lose a large 
fraction of rotation velocity, but only for one galaxy (GAL6) 
we do not observe any significative rotation along the major 
axis. GAL8 shows a well defined symmetric rotation pattern 
although its peak velocities is only ~ 7 km s _1 . GAL6 and 
GAL8 do not show evidences for structures or disk features 
and are spheroidal systems with a c/a ratio larger (Table 1) 
than the other remnants and surface brightness profiles well 
fitted by exponential laws. Moreover, they are quite small 
(M s 1.5 x 10 9 Mq) baryonic dominated galaxies presently 
orbiting at ~ 100 kpc from the cluster center. 

Is the lack of rotation observed in GAL6 due to a pro- 
jection effect, i.e. we are looking at the face on projection 
of a rotating system, or is it physically significative? Fig. 
16 illustrates the kinematic profiles of GAL6 and GAL8 for 
different inclinations of the remnant: from the top to the 
bottom i — 0°, which corresponds to a face-on projection, 
i = 45° and i — 90° (edge-on). The last panel on the bottom 
represents an end-on galaxy, with inclination i = 90° and the 
three dimensional major axis aligned with the line of sight. 



GAL6 shows traces of rotation only for i = 0, with a maxi- 
mum rotational velocity v ma x = 3.7kms _1 , while the other 
projections have not well defined rotation curves and val- 
ues of Vmax^j 2.7 km s -1 . In the case of non-rotating galaxies 
we estimated an upper limit for v max by differentiating the 
average vel ocities on each side of the major axis and di vid- 
ing by two iGeha. Guhathakurta fc van der Marel20"02Tl . In 
the calculation of v max for i = 0° and i = 90° we excluded 
points with velocity dispersion larger than 80kms _1 , repre- 
senting stars close to the tidal radius of the remnant, which 
projected motion respect to the center of mass is due to 
stripping processes and not to rotation. On the other hand 
GAL8 has always a well defined rotation pattern and there 
is not a preferential direction with zero rotation. Among the 
other galaxies characterised by low rotation velocities (see 
Fig. 15), GAL2, GAL9 and GAL10 are more massive sys- 
tems if compared with the two spheroids GAL6 and GAL8, 
but have still high c/a values, while the low line of sight 
velocity of GAL12 is not related to a morphological trans- 
formation but is simply due to the fact that it is a spiral 
galaxy viewed face on (Fig. 14). 

On the basis of this limited sample we expect non- 
rotating and slowly rotating dwarfs to have internal and 
environmental properties substantially different from the ro- 
tating ones. In particular the former seem to be the product 
of a violent morphological transformation accelerated by the 
proximity of the cluster center. This result is apparently in 
contradiction with Geha et al. 2003 who finds no relation- 
ship between the amount of rotation and the morphology or 
the orbital properties of the galaxies in their sample. 

From Fig. 15 it appears evident that although the galax- 
ies lose significant amounts of mass, the velocity dispersion 
can rise up to 50kms _1 in the central 1-2 kpc. This is due 
to the bar formation and the transformation of circular to 
more radial orbits. The velocity dispersion profiles are typ- 
ically decreasing at larger radii but they rise again in those 
galaxies for which the slit semi-length approaches or exceeds 
the tidal radius. 

As the ratio between rotational velocity and velocity 
dispersion decreases, the galaxy loses rotational support and 
the flattening becomes dominated by anisotropic pressure. 
The degree of rotational support can be expressed in terms 
of the anisotropy parameter (v/a)* llBinnevlll978t) . defined 
as 



where v is the rotational velocity, a and e the mean projected 
velocity dispersion and ellipticity. The theoretical prediction 
for a rotationally flattened spheroidal is given by (v/a)* ~ 1 
jBinnev fc Tremaindll"987h and is represented in Fig. 20 by 
a continuous solid line. Galaxies well below this curve have 
significant velocity anisotropy. 

In order to compare our results with observations, we 
plotted in Fig. 17 the ratio of the maximum rotational veloc- 
ity Vmax to the mean velocity dispersion a within a radius 
of 6 kpc versus the mean isophotal ellipticity calculated in 
the same radial range. In the same plot we show the obser- 
vational results for most of the dwarf galaxies in Virgo with 
published kinematics. We a lso plot the data for dw arf galax- 
ies in the Fornax cluster bv lDe Riicke et all l)2003|) . Some of 
the galaxies in Fig. 17 (VCC 543, VCC 917, VCC 1036, VCC 
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Figure 17. The ratio of the maximum rotational velocity v max 
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1261 and VCC 1308) were observed by both Geha et al. 2003 
and van Zee et. al 2003, using the Keck II and Palomar tele- 
scope. The observations of van Zee et. al 2003, which cover a 
larger radial extent due to the aperture size, gives to signif- 
icantly higher values of v max respect Geha et al. 2003. This 
appears to be due to the fact that the smaller slit used by 
the latter authors only samples the inner rising part of the 
rotation curve and does not include the maximum value. For 
example VCC 917 and VCC 1308, classified as non-rotating 
galaxies (both with v ma x = 0.4 km s _1 ) by Geha et al. 2003, 
have Vmax > 10 km s _1 according to van Zee et al. 2003, 
while the rotation curves of the slowly rotating dwarfs VCC 
543 and VCC 1036 (v max = 13 and 14 km s" 1 in Geha et 
al. 2003) peak beyond 15" with v ma x > 40 km s . Most of 
the simulated remnants have a similar range of ellipticities 
and Vmax/ (J as the observed dwarfs, even if it is clear that 
we have some difficulties in reproducing galaxies with a de- 
gree of rotational support close to zero. On the other hand 
it is possible that the low v ma x/a ratio of some galaxies is 
due to the small radial range covered by the observations. 
In order to understand how much the slit length influences 
the results, in Fig. 18 we plotted the v/a ratio of the pro- 
jected remnants as a function of the radial distance along 
the slit. In most of the cases v/a increases with radius (with 
the remarkable exception of GAL1 whose rotational support 
drops beyond r = 2.5 kpc). Using a 4 kpc slit we would find 
10 galaxies below v m ax/a = 1, 6 of which with » m „/a < 1 
and 3 well below v m ax/a = 0.1, which is the upper limit 
assumed by Geha et al. 2003 for non-rotating dwarfs. In- 
creasing the slit length up to 8 kpc, the number of galaxies 
with a low degree of rotational support is still quit high (7 
galaxies with v m ax/a < 1 and 4 with v m ax/a < 0.5), but 
only GAL6 has still a v ma x/a ~ 0, as in the case of the 
adopted slit length of 12 kpc. 




Figure 18. v/a ratio as a function of the radial distance along 
the major axis for the 13 projected remnants of Fig. 15. 



6 CONCLUSIONS 

We follow the evolution of late type disk galaxies as they 
evolve within a hierarchically forming galaxy cluster using 
N-body simulations. We use a single galaxy model to explore 
the effects of cluster-centric position and orbit. Our resolu- 
tion is such that we can compare with Keck and Palomar 
kinematics and VLT imaging data. We find that: 

(i) The cluster environment can be responsible for trig- 
gering bar and buckling instabilities. 

(ii) Tidal shocks with dark matter substructures and 
the mean cluster potential are both important at driving 
the evolution from disks to spheroidal structures. However, 
galaxies within the inner 100 kpc suffer the most complete 
transformation. 

(iii) None of the disks are completely destroyed there- 
fore we do not expect to cre ate UCD galaxies b y progenitor 
disks with cuspy dark halos llBekki et alJl2003l) . 

(iv) More tidal heating leads to more mass loss which 
leads to rounder galaxies. 

(v) The kinematics (v/a) of the stellar remnants is sim- 
ilar to observations of Virgo and Fornax dE/dSph galaxies. 
Only the most severely disrupted galaxies which lose most 
of their stellar mass become completely pressure supported. 
Most of the remnants retain a significant amount of rota- 
tional motion. 
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(vi) Unsharp masked images of the galaxies show sim- 
ilar hidden bar and spiral features as some cluster galax- 
ies. These gradually dissapear as the galaxies suffer more 
heating. Our results could be rescaled in time and mass 
to make statements about possible disk progenitors enter- 
ing group or even galactic environments. The next step is 
to self-consistently follow the evolution of galaxies as they 
form and evolve in different environments. 

Our study has focussed on the evolution of a late type 
disk similar to the Large Magellanic Cloud. The baryon frac- 
tion in the disk is relatively high and adiabatic contraction 
has enabled the galaxy to be more robust. Disks such as 
M32, have a baryon fraction about five times lower and thus 
will be more susceptible to tidal perturbations. Exploring 
a large parameter space is beyond the scope of this paper, 
however we note that if we followed the evolution of fainter, 
baryon poor disks we would expect a more complete trans- 
formation to dE/dSph and a lower fraction of galaxies with 
hidden features such as bars and disks. 



7 ACKNOWLEDGEMENTS 

We would like to thank Jiirg Diemand for providing 
the cosmological cluster simulation. We are grateful to 
Fabio Barazza, Maria Geha, Frank van den Bosch, Ste- 
lios Kazantzidis and Andrea Maccio for useful discussions. 
The numerical simulations were performed on the Zbox 
(http:/www.theorie.physik.unizh.ch/stadel) supercomputer 
at the University of Zurich. CM and RP are supported by 
the Swiss National Science Foundation. VPD is supported 
by a Brooks Fellowship at the University of Washington. 



REFERENCES 

Barazza, F. D., Binggeli, B., Jerjen, H. 2002, A&A, 391, 
823 

Barazza, F. D., Binggeli, B., Jerjen, H. 2003, A&Ap, 407, 
121 

Bckki, K., Couch, W. J., Drinkwater, M. J., Shioya, Y. 

2003, MNRAS, 344, 399 
Bender, R., Nieto, J. L. 1990, A&A, 239, 97 
Binney, J. 1978, MNRAS, 183, 501 

Binney, J., Trcmainc, S. 1987, Princeton, NJ, Princeton 

University Press, 1987 
Conselice, C. J., Gallagher, J. S., Wyse, R. F. G. 2001, 

ApJ, 559, 791 
Courteau, S. 1997, AJ, 114, 2402 

De Rijcke, S., Dejonghe, H., Zeilinger, W. W., Hau, G. K. 
T. 2001, ApJ, 559, L21 

De Rijcke, S., Dejonghe, H., Zeilinger, W. W., Hau, G. K. 
T. 2003, A&A, 400, 119 

Diemand, J., Moore, B., Stadel, J. 2004a, MNRAS submit- 
ted (astroph 0402267) 

Diemand, J., Moore, B., Stadel, J. 2004b, MNRAS, tmp. 
152D 

Dressier, A., Oemler, A. J., Butcher, H. R., Gunn, J. E. 

1994, ApJ, 430, 107 
Ferguson, H. C, Sandage, A. 1991, AJ, 101, 765 
Freedman, W. L., et al. 2001, ApJ, 553, 47 
Gavazzi, G. et al, astroph 0410228 Accepted to A&A 



Geha, M., Guhathakurta, P., van der Marel, R. P. 2002, 
AJ, 124, 3073 

Geha, M., Guhathakurta, P., van der Marel, R. P. 2003, 

AJ, 126, 1794 
Gnedin, O.Y. 2003, ApJ, 589, 752 

Graham, A. W., Jerjen, H., Guzman, R. 2003, AJ, 126, 
1787 

Hernquist, L. 1993, ApJ, 86, 389 

Jerjen, H., Kalnajs, A. & Binggeli, B. 2000, A&A, 358, 845 
Kormendy, J., Freeman, K. C. 2004, IAU Symposium, 220, 
377 

Mayer, L., Governato, F., Colpi, M., Moore, B., Quinn, T., 
Wadsley, J., Stadel, J., Lake, G. 2001, ApJ, 559, 754 

Mayer, L., Moore, F., Quinn. T., Governato, F., Stadel, J., 
2002, MNRAS, 336, 119 

Mo, H. J., Mao, S., White, S. D. M. 1998, MNRAS, 295, 
319 

Moore, B., Katz, N., Lake, G., Dressier, A., Oemler, A. 

1996, Nature, 379, 613 
Moore, B., Lake, G. & Katz, N. 1998, ApJ, 495, 139 
Moore, B., Lake, G., Quinn, T. & Stadel, J. 1999, MNRAS, 

304, 465 

Moore, B., Diemand, J., Stadel, J. 2004, IAU Colloquium: 
Outskirts of Galaxy Clusters: intense life in the suburbs 
(astroph 0406615) 

Navarro J. F., Frenk C. S., White S. D. M. 1996, ApJ, 462, 
563 

Navarro J. F., Frenk C. S., White S. D. M. 1997, ApJ, 490, 
493 

O'Neil, K., Bothun, G. D., Impey, C. D. 1999, AJ, 118, 
1618 

Pedraz, S., Gorgas, J., Cardiel, N., Sanchez-Blazquez, P., 

Guzman, R. 2002, MNRAS, 332, L59 
Persic, M., Salucci, P. 1997, ASP Conf. Ser. 117: Dark and 

Visible Matter in Galaxies and Cosmological Implications, 
Simien, F., Prugniel, P. 2002, A&A, 384, 371 
Sparke, L. S., Sellwood, J. A. 1987, MNRAS, 225, 653 
Spergel, D. N. et al., 2003, ApJSS, 148, 175 
Springel, V., White, S. D. M. 1999, MNRAS, 307, 162 
Stadel, J. 2001, Ph.D. Thesis 
Toomrc, A. 1964, ApJ, 139, 1217 

Toomrc, A. 1981, Structure and Evolution of Normal 
Galaxies, S. M. Fall & D. Lynden-Bell, Cambridge Univ. 
Press, 111 

van Zee, L., Skillman, E. D., & Haynes, M. P. 2004, AJ, 
128, 121 

van Zee, L., Barton, E. J., Skillman, E. D., astroph 0409346 
Accepted to AJ. 



